Single spin- and chiral-glass transition in vector spin glasses in three-dimensions 
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Results of Monte Carlo simulations of XY and Heisenberg spin glass models in three dimensions 
are presented. A finite size scaling analysis of the correlation length of the spins and chiralities 
of both models shows that there is a single, finite-temperature transition at which both spins and 
chiralities order. 
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There is now general consensus that a spin glass tran- 
sition occurs in three dimensional Ising spin glasses with 
short range interactions at finite temperature The 
most convincing work is that of Ballesteros et al. jl| who 
performed a finite-size scaling analysis of the correlation 
length, ^i, in samples of different sizes L. Data for the 
dimensionless ratio ^l/L is found to intersect cleanly at 
T = TsGy as expected at a second order phase transition. 

The situation is much less clear, however, for vec- 
tor spin glasses. Early work on XY 0, 01 and Heisen- 
berg 0, 0] models indicated a zero temperature tran- 
sition, or possibly a transition at a very low but non- 
zero temperature. However, following earlier work of Vil- 
lain , which emphasized the role of "chiralities" (Ising- 
like variables which describe the handedness of the non- 
collinear spin structures), Kawamura and Tanemura 
proposed a chirality transition at T = Tea (> 0), even 
though the spin glass transition temperature, Tgc, is as- 
sumed to be zero. This scenario requires that spins and 
chiralities decouple at long length scales. Kawamura and 
collaborators have given numerical evidence for this sce- 
nario both for XY Q and Heisenberg H,0] models. 

However, the absence of a spin glass transition in 
vector spin glass models has been challenged. For the 
XY case, Maucourt and Grempel 0| and subsequently 
Akino and Kosterlitz found evidence for a possible 
finite Tsg from zero temperature domain wall calcula- 
tions. Furthermore, by studying the dynamics of the XY 
spin glass in the phase representation, Granato |3| found 
that the "current-voltage" characteristics exhibited scal- 
ing behavior which he interpreted as a transition in the 
spins as well as the chiralities. For the Heisenberg model, 
Matsubara et al. Q| , and Nakamura and Endoh 
have argued that the spins and chiralities order at the 
same low but finite temperature. 

Since the most successful approach to demonstrate a 
finite temperature transition in the Ising case has been 
the scaling of the correlation length il| it seems useful to 
perform a similar analysis for vector spin glass models. 
Furthermore, one can calculate the correlation length for 
both the spins and chirality, and so perform the same 
analysis for both types of ordering. Here, we present re- 
sults of these calculations for the XY and Heisenberg 
models. For both models, we find a single transition for 



both spins and chiralities at low but finite temperature. 

We take the standard Edwards- Anderson spin glass 
model 
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where the Si are n-component vectors of unit length at 
the sites of a simple cubic lattice, and the Jy are near- 
est neighbor interactions with zero mean and standard 
deviation unity. We consider both the XY model (n—2), 
and the Heisenberg model (n = 3). Periodic boundary 
conditions are applied on lattices with N = spins. 

The spin glass order parameter generalized to wave 
vector k, (/^''(k), is defined to be 



(2) 



where /i and i' are spin components, and "(1)" and "(2)" 
denote two identical copies of the system with the same 
interactions. From this we determine the wave vector 
dependent spin glass susceptibility xsgC^) by 



XsG(k)=iVV[(k'^'^(k 



(3) 



where (• • •) denotes a thermal average and [• • -Jav de- 
notes an average over disorder. The spin glass correlation 
length is then determined P,^^, from 



Xsg(O) 
2sin(fcmin/2) ^ 

V XSG (krniii) 
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where k,„i„ = (27r/L)( 1,0,0). 

For the XY model the chirality of a square is |3 
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sgn(J;,„) sin(6'; - 6m), 
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where 9i is the angle characterizing the direction of spin 
S/, and the sum is over the four bonds around the ele- 
mentary plaquette perpendicular to the ^ axis and whose 
"bottom left" corner is site i. The chiral glass suscepti- 
bility is then given by 



(6) 
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FIG. 1: A plot of the spin (^l) and chiral (^jj ^ and fi;^) 
correlation lengths for different sizes and temperatures for 
the XY spin glass. The solid lines connect data for the 
dotted lines connect ^^^i^, and the dashed lines connect ^| j^. 



FIG. 2: Plot of the spin glass correlation length divided 
by L for the XY spin glass. The data intersects at T ~ 
0.34, implying that there is a spin glass transition at this 
temperature. The inset shows a scaling plot according to 
Eq. lO with TsG = 0.33 and v = 1.2. 



where the chiral overlap q^(k) is given by 

i 

We define the chiral correlation lengths by 
^ I ( XccjO) A'^' 

^^'^ 2sin(fc,ni„/2) Vx&G(kmin) / ' ^ ' 

in which xcgO^ — 0) is independent of fi. Note that 
L will, in general, be different for fl along kmin (the x 

direction) and perpendicular to k. We denote these two 

lengths by ^ and respectively. 

For the Heisenberg spin glass, Kawamura defines 

the local chirality in terms of three spins on a line as 

follows: 

l^i = ^i+fi ' X Sj_^. (9) 

The chiral glass susceptibilities and correlation lengths 
are then given in terms of the /tf by Eqs. 0-®, as for 
the XY model. 

We use parallel tempering 0, 0| Monte Carlo to go 
down to the low temperatures that are needed, and study 
sizes from L = 4 to 12. To test fo r eq uilibration we 
require that the following relation [20j, 

[qi-qsU^-TpU, (10) 



valid for a Gaussian bond distribution, is satisfied. Here 
U is the energy per spin, qi = {l/Nb)J2{ij){^i ' is 
the "link overlap", q, = (IM) ((s/ • Sj)^) where 

Nh — {z/2)N is the number of nearest neighbor bonds, 
and z 6 here) is the lattice coordination number. 
We averaged over 1000 samples, except for the follow- 
ing cases: XY, L = 12, 601 samples, Heisenberg: L = 8, 
436 samples, and L = 12, 331 samples. The number of 
sweeps that each set of spins performed varied from 6000 
for the small sizes to 300000 for L = 12. 

Since /L is dimensionless it has the finite size scaling 
form 

^=x(l1/^(T-Tsg)), (11) 

where v is the correlation length exponent. Note that 
there is no power of L multiplying the scaling function 
X, as there would be for a quantity with dimensions. 
There are analogous expressions for the chiral correlation 
lengths. From Eq. (|ll|l it follows that the data for £,l/L 
different sizes come together at T = TsG- In addition, 
they are also expected to splay out again on the low-T 
side if there is spin glass order below Tsg- 

Next we discuss the results, starting with the XY spin 
glass. Data for the various correlation lengths are shown 
in Fig. n One sees that the chiral correlation lengths 
are smaller than the spin glass correlation length at the 
higher temperatures, but increase faster on lowering T, 
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FIG. 3; The main figure shows data for the perpendicular 
chiral correlation length ^^f^ divided by L, for sizes 4 < L < 
12 for the XY spin glass. There are intersections at about the 
same temperature as that found for the spins in Fig. |21 but 
evidently with some corrections to scaling. The inset shows 
analogous data for the parallel chiral correlation length, ^jj ^ . 

such that, for a given L, all the lengths become compa- 
rable at the lowest temperatures simulated. Furthermore 
the two chiral correlation lengths (parallel and perpen- 
dicular to k) become indistinguishable at lower T and 
larger sizes, as one would expect. 

The data for S^l/L, shown in Fig. [2 intersects at a 
well defined temperature ~ 0.34 and splay out at lower 
temperatures, which, according to Eq. implies a 

transition at this temperature. We find 

TsG = 0.34 ± 0.02 XY spin glass. (12) 

Note the L = 12 data intersects at somewhat lower T, 
implying that corrections to finite size scaling may still 
be significant for this range of sizes. Fig. El provides com- 
pelling evidence, in our view, that there is finite spin 
glass transition temperature in a three-dimensional XY 
spin glass, in contrast to the claim in most of the liter- 
ature. The inset to Fig. [21 shows the data collapses well 
according to Eq. pi|l with i/ = 1.2 ± 0.2. In this paper, 
error bars do not include systematic effects which are 
hard to estimate. Given the lower intersection point of 
the L = 12 data, it is possible that Tc could be lower than 
that estimated here, in which case the value of v would 
be increased perhaps to the Ising valueflj 2.15 ± 0.15. 

Data for the chiral correlation lengths for the XY 
model is shown in Fig. Though the intersections are 
not quite as clean, they occur at about the same tempera- 
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FIG. 4: Data for the spin glass correlation length ^l, divided 
by L for the Heisenberg spin glass. The intersections imply 
that TsG — 0.16. The inset shows a scaling plot according to 
Eq. dJ with TsG = 0.16 and v = 1.1. 

ture as found above for the spin glass correlation length. 
For both spin and chiral correlations, they occur at a 
slightly lower T for the L ^ 12 data. Altogether, the ev- 
idence is good that there is a finite chiral glass transition 
at (or very close to) TsG- Collapsing the data we find 
the chiral correlation length exponent is Vc = 1.3 ± 0.3, 
which is compatible with our estimate for the spin cor- 
relation exponent ly. Note that if the spins order then 
the chiralities must also order, assuming a non-coUinear 
state, and so Tcg > Tsg- 

We are not aware of any estimates of transition tem- 
peratures for the XY spin glass with Gaussian couplings, 
though for the ±J model Kawamura and Li [7| find 
Tcg — 0.39 ± 0.03, somewhat higher than ours. Since, 
for the Ising spin glass, Tc is somewhat higher for the 
± J model than the Gaussian model, our estimate of Tc 
is probably compatible with Kawamura and Li's. How- 
ever, we emphasize that, in contrast to them, we find 
simultaneous ordering of the spins and chiralities. 

Next, we go on to our results for the Heisenberg spin 
glass. As for the XY model, the spin glass correlation 
length is larger at higher temperatures but the chiral 
correlation length grows faster and is comparable to the 
spin correlation length at the lowest temperatures. Fig- 
ure 21 shows data for S,l/L, which intersect at a common 
temperature indicating a finite spin glass transition tem- 
perature Tsg which we estimate to be 

Tsg = 0.16 ± 0.02 Heisenberg spin glass. (13) 
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Why has this simple picture of a single finite tempera- 
ture transition in vector spin glasses in three dimensions 
not been generally accepted before? One reason is that 
TsG is very low compared with the mean field value T^q^ 
(~ f.22 for XY and 0.82 for Heisenberg). Until the ad- 
vent of parallel tempering 0| it was difficult to reach the 
actual TsG in simulations. Furthermore, the comm only 
used Binder ratio does not seem to be very useful |21| 
for vector spins, and T = domain wall calculations 
are plagued by uncertainties over the optimal choice of 
boundary conditions We argue that, as for the Ising 
spin glass , finite size scaling of the correlation length 
is the optimal technique to use. 

We would like to thank Helmut Katzgraber for helpful 
discussions, and both he and H. Kawamura for comments 
on a preliminary version of the manuscript. We acknowl- 
edge support from the National Science Foundation un- 
der grant DMR 0086287. 



FIG. 5; The main figure shows data for the parallel chiral 
correlation length divided by L for the Heisenberg spin 
glass. There is an intersection at T ~ 0.14, close to that found 
for the spin glass correlation length in Fig. |1] The inset 
shows analogous data for the perpendicular chiral correlation 
length. The data intersect but not as cleanly as for ^^f^ or ^l- 



The inset to Fig. ^ shows the data collapses quite well 
according to Eq. (^TJ with v = 1.1 ± 0.2. 

Figure. [SI shows that the data for the chiral correlation 
lengths indicate a transition at about the same value, 
with the intersections being cleaner for the parallel than 
for the perpendicular correlation length. Our estimate for 
the chiral correlation length exponent is Vc — 1.3 ± 0.3. 
As for the XY spin glass, Tc may be somewhat lower than 
that found here, which would lead to a larger value of v. 

Our value for the transition temperature agrees well 
with values of Tcg given by Kawamura "il, 0.157±0.01, 
and Hukushima and Kawamura 9], 0.160±0.005, though, 
unlike those authors, we claim that the spins, as well as 
the chiralities, order at this temperature. For the ±J 
model, Endoh et al. |Q find a spin glass transition at 
TsG 0.19 ± 0.02 while Nakamura and Endoh ^ find 
both chiral and spin glass ordering for T ~ 0.21. 

To conclude, by analyzing data for the spin and chiral 
correlation lengths, we have argued that there is a sin- 
gle phase transition, at which both spins and chiralities 
order, in the XY and Heisenberg spin glasses in three 
dimensions. In our view, the evidence for a spin glass 
transition is at least as strong as that for a chiral glass 
transition. Spin-chirality decoupling does not seem to 
occur. The present work used quite modest work station 
facilities, so it would be feasible to extend these results 
to larger sizes by a major computational effort. 
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